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Epidemic Incidence in Correlated Complex Networks 
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We introduce a numerical method to solve epidemic models on the underlying topology of complex 
networks. The approach exploits the mean-field like rate equations describing the system and allows 
to work with very large system sizes, where Monte Carlo simulations are useless due to memory 
needs. We then study the SIR epidemiological model on assortative networks, providing numerical 
evidence of the absence of epidemic thresholds. Besides, the time profiles of the populations are 
analyzed. Finally, we stress that the present method would allow to solve arbitrary epidemic-like 
models provided that they can be described by mean-field rate equations. 
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A few years ago, Watts and Strogatz [l| introduced a 
model able to produce networks with properties of both 
regular lattices and random graphs with small diameter. 
Their model soon led to a burst of activity in the field 
0, El , further spurred by Barabasi and collaborators who 
found that many seemingly diverse systems share sev- 
eral topological properties such as a power law behavior 
in their connectivity distributions when represented as 
networks These complex networks are formed by a 
set of many elements (or nodes) that are linked together 
through edges (or links) if they interact directly. Empir- 
ical evidence supports that in notable networks, such as 
metabolic or communication webs, the probability P(k) 
that any node has k links to other nodes is distributed 
according ly to a power law P(k) ~ k~~< 0, H 0, 

with 

7 < 3 in most cases. 

Networks of this type, called scale-free (SF) networks, 
show a noticeable property: the heterogeneity of the con- 
nectivity distribution can not be neglected. One of the 
fundamental results derived from this property is that 
the threshold characterizing the percolation transition 
or an epidemic outbreak is vanishing in the thermody- 
namic limit 0, 0, 0|. On the other hand, many 
real- world networks are also characterized by degree cor- 
relations [eH Il3| , which make it necessary to reconsider 
the same problems but taking into account the condi- 
tional probability P{k'\k) that a link emanating from 
a node of connectivity (or degree) k leads to a node 
of connectivity k' . Very recently, it has been shown 
analytically that the presence of nontrivial correlations 
does not change the main conclusions drawn for uncor- 
rected graphs, namely, the absence of percolation and 
epidemic thresholds [Til under very general condi- 
tions. However, neither complete exact solutions nor nu- 



merical studies of the system's dynamics have been re- 
ported. From a theoretical side, it seems an unsurmount- 
able task to work out such study. In contrast, although 
feasible in principle, numerical simulations using Monte 
Carlo techniques would be quite arduous since one should 
first generate a network with the proper correlations and 
then perform expensive Monte Carlo simulations for large 
system sizes, where memory requirements could seriously 
limit the size of the networks under study. 

In this paper, we present an efficient numerical method 
that allows the study of epidemic-like models on the un- 
derlying topology of complex networks with arbitrary 
connectivity distributions and degree-degree correlations. 
The method solves the mean-field rate equations describ- 
ing the dynamics of the system, where the topological 
properties are accounted for. Specifically, we fully an- 
alyze the Susccptiblc-Infectcd-Rcmoved (SIR) model in 
networks with assortative correlations, where nodes tend 
to be linked with their connectivity peers as for the 
case of social networks. Our results are compared with 
those obtained for uncorrelated networks and confirm 
the previous findings about the absence of any epidemic 
threshold. Besides, we report on the time evolution of the 
individual populations which allow us to draw interesting 
conclusions when confronted to uncorrelated networks. 

The SIR model [HG3 considers that individuals are 
classified in three classes according to their state: suscep- 
tible, infected and removed. The epidemic is propagated 
by contacts between infected and susceptible individuals 
at a rate A. Thus, once an individual gets infected and re- 
covers he can not catch the epidemic again. Moreover, in 
networks with power-law distributed connectivities one 
also has to consider the presence of nodes with differ- 
ent connectivity k within each category. We consider 



the time evolution of the magnitudes p k (t), s k (t), and 
rfe(t), which are the density of infected, susceptible, and 
removed nodes of connectivity k at time t, respectively, 
with the normalization condition p k (t)+s k (t)+r k (t) = 1. 
Global quantities can be obtained by averaging over the 
connectivity classes. In this way, the fraction of in- 
fected individuals (or the epidemic incidence) is given 

ttfK*)=EfcW-*(t). 

The mean-field rate equations for the evolution of these 
densities satisfy the following set of coupled differential 
equations |17|: 



ds k (t) 
dt 

dpk{t) 
dt 

dr k {t) 
dt 



-Xks k (t)J2 p {k'\k)Pk', (1) 
k> 

-p k (t)+\ks k (t)J2 p ik'\k)pk>, (2) 

k' 

Pk(t). (3) 



with the initial conditions s(0) = (N — l)/N, p(0) = 1/N 
and r(0) = 0. In the above equations, it is consid- 
ered that infected individuals recover with unitary rate. 
Moreover, the creation term in Eq. @ is proportional 
to the density of susceptible individuals, s k (t), times the 
spreading rate A, the number of emanating links k and 
the probability that any neighboring node is infected. 
This probability is given by the average over all degrees 
of the probability P{k'\k)p k i that a link emanated from a 
node with connectivity k points to an infected node with 
degree k' [UGH- 

The numerical approach introduced here is based on a 
different interpretation of the mean-field rate equations 
(|1I3|I and is implemented as follows. Indeed, these equa- 
tions describe a process in which individuals are decaying 
from one state to another. Hence, one can speak in terms 
of transition probabilities from the class of susceptible 
onto infected and finally onto removed individuals, i.e., 

Sk — ► Pk — ► T k . From Eqs. 1113(1 we get the transition 
probabilities at each time step t: 

W* p {t) = XkNP(k) Sk (t)J2P(k'\k)pk'(t), (4) 



W k pr {t) = NP{k)p k {t), 



(5) 



where all the topological information is contained. The 
mean time interval, r, for one transition to occur after 
i — \ decays is, 



1 



w sp {t) + w pr {ty 



(6) 



with W sp (t) = £ fc W*(t), W pr (t) = EkKrV and 
t = X T j- That is, at any instant t of the decay- 
ing process, r represents the mean time for the next in- 
dividual decay. Once the transition probabilities 14151) 
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FIG. 1: Epidemic incidence in scale- free networks with 7 = 3 
and TV = 10 5 with (full squares) and without correlations (full 
circles) as a function of the spreading rate A. The assortative 
correlations are given by Eq. Q with a = 0.1. Fhe inset is 
a fit to the analytic relation R ~ e _2 /( fc >\ See the text for 
details. 



are calculated, we stochastically decide what transition 
takes place. Nodes are divided in three classes according 
to their state and, within each of these classes, they are 
also characterized by their connectivity k. Hence, the 
identification of both what transition occurs and which 
class k is affected after one r is done by deciding that 
the probabilities that precisely a node with connectivity 
k changes its state at time t are given by 



K^) = W* (i)r, 



II* (t) = W* (t)r, (7) 



materializing the choice by generating a random number 
between and 1. In other words, the individual decays 
within each connectivity class k proceeds by chance with 
the probabilities dictated by Eq. until the end of the 
spreading process, i.e., when the condition p k — is veri- 
fied Vfc. It is worth mentioning that the present stochas- 
tic approach has been successfully applied to the study 
of time-dependent models of fracture [20I l2l|. Besides, 
we remark that the main advantage of the method lies in 
the fact that we don't have to know explicitly the con- 
nectivity matrix, but only a functional form for P(k'\k). 
Therefore, in contrast to Monte Carlo simulations, we 
do not generate any network. Instead, we generate a se- 
quence of integers distributed according to a power law 
P(k) ~ k^ 1 for the node's connectivities. 

In order to compare the SIR dynamics in uncorrelated 
and correlated networks, we set henceforth 7 = 3. For 
uncorrelated networks, the two-point degree correlation 
function P(k'\k) is of the form 



P(k'\k) = qy 



k'P(k') 



(8) 
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FIG. 2: Likelihood of an epidemic outbreak as a function of 
the connectivity of the initial seed. The value of A has been 
set to 0.15. Note that the introduction of correlations does 
not change the probability of having an outbreak. The figure 
is in log-linear scales for clarity. The network's parameters 
are as of Fig. Q 

Furthermore, consider the case in which the degree cor- 
relations can be decomposed into two components 

P{k'\k) = (l-a)q k , +aS kk ,, (9) 

with < a < 1. Varying the parameter a one inter- 
polates between the uncorrelated graphs (a = 0) and a 
graph with positive degree correlations 0, • 

We have studied the effects introduced by assortative 
correlations given by Eq. J^J in the spreading of a dis- 
ease. Social networks are the capital example of assor- 
tative networks and, although positively correlated, do 
not have high correlation coefficients [lfij |. hence, we have 
used small values of a in our simulations. Figure ^ shows 
the epidemic incidence, R — r(t oc ) = ^2 k P{k)rk{too), 
averaged over many realizations, where too is the lifetime 
of an epidemic outbreak. The results indicate that the 
epidemic incidence in networks with assortative correla- 
tions is smaller that for uncorrelated ones with the same 
degree distribution. However, the epidemic threshold in 
the thermodynamic limit is not modified by the presence 
of positive correlations and is vanishing when N — > oo. 
The inset in Fig. ^ is a plot of R as a function of 1/A 
for an uncorrelated network with (k) — 6. The behavior 
nicely fits the theoretical prediction R ~ e _2 /^ A [l7j . 
being the slope of the straight line equal to 0.36(1). This 
constitutes a further test of the validity of the approach 
presented here. As to the correlated networks, they sat- 
isfy the same functional form R ~ e~ c l x but approaches 
zero more slowly. Thus, the numerical findings confirm 
some theoretical arguments recently pointed out about 
the absence of any epidemic threshold in SF correlated 



FIG. 3: Epidemic incidence as a function of the connectivity 
of the initial seed for A = 0.15. The dependency on the initial 
seed is stronger for uncorrelated networks. The network's 
parameters are as of Fig. Q See the text for further details. 



networks with 7 < 3 |l4. Il5|. Moreover, we note that in 
finite size networks with assortative correlations the ef- 
fective threshold is larger, suggesting that these networks 
are more robust than uncorrelated ones. 

The high heterogeneity of SF networks also causes that 
the relative incidence of an outbreak stron gly d epends on 
the connectivity of the first infected nodes 17] . First, we 
explored the likelihood of an epidemic outbreak, P out as 
a function of the connectivity, ki, of the initial infected 
node. This probability is obtained by dividing the num- 
ber of times an epidemic developed by the total attempts 
made. The results are summarized in Fig. [5J As can 
be observed, no matter whether or not correlations are 
present, the probability of having an epidemic outbreak 
is the same for both kind of networks. This behavior 
changes when coarse-graining the results of Fig. 2] The 
results shown in Fig. have been obtained by simulating 
the SIR dynamics when the initial infected node has a 
connectivity ki and recording the epidemic incidence in 
each case. Now, the number of individuals within the 
removed class at the end of the spreading process sig- 
nificantly depends on the correlations. Note that this is 
not a consequence of the smaller value of the epidemic 
incidence for assortative correlations. In fact, the first 
part of the curves, up to intermediate values of ki, fol- 
lows a logarithmic dependency on the connectivity, but 
the number of removed individuals for uncorrelated net- 
works grows more faster than for correlated ones. This 
behavior can be intuitively understood since in assorta- 
tive networks, nodes are mainly connected with their con- 
nectivity peers (even for small values of a), whereas in 
uncorrelated networks it is always possible that a poorly 
connected node transmits the disease to a hub. In other 
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FIG. 4: Time evolution of the epidemic incidence (main fig- 
ure) and of the density of infected individuals (inset) when the 
initial infected individual has the largest connectivity k ma x for 
both uncorrelated and correlated networks. A is equal to 0.15. 
The network's parameters are as of Fig0 



words, the random mixing leads to a higher R and favors 
the propagation of the outbreak. 

The method introduced here also allows the explo- 
ration of the time dependency of the quantities character- 
izing the epidemic spreading. Note that the time profiles, 
in contrast to MC simulations, are in units of A~ x and 
thus they can be used as a real quantitative measure. Fig- 
ure 01 plots the time evolution of the fraction of removed 
individuals and that of the infected nodes for the two 
networks under study when the initial infected individ- 
ual has the largest connectivity k max . Here, as before, by 
simple inspection one would not be able to distinguish the 
behavior of these magnitudes in both networks as they 
show the same functional dependency on time. However, 
a more careful look at the plateau of R(Xt) reveals that 
the epidemic lifetime is longer for assortative correlations 
than for uncorrelated networks. This is a direct conse- 
quence of the correlations. Starting from the very hub 
of the network, as time goes on, the epidemic is prop- 
agated, on average, from highly connected nodes down 
to less connected individuals such that when the end of 
the process is approaching, the mean time for individual 
decays is longer and longer leading to an effective de- 
celeration of the epidemic spreading process. This is not 
anymore the case when random mixing rules out any cor- 
related spreading. A further evidence of this mechanism 
is provided in the inset where it can be clearly noted that 
the two density profiles crosses well before the final death 
of the disease. 

In summary, we have introduced an efficient numerical 
method that allows to explore the spreading of epidemic 
diseases in complex correlated networks without gener- 



ating explicitly the network with the proper correlations. 
We have studied the SIR epidemiological model in as- 
sortative networks and found that its qualitative behav- 
ior is the same as for uncorrelated networks. However, 
from a more practical perspective, there are some im- 
portant quantitative differences that deserve to be con- 
sidered more carefully and makes the inclusion of cor- 
relations in nowadays studies meaningful. In particular, 
while the likelihood of an epidemic outbreak is not mod- 
ified when taking into account positive correlations, the 
epidemic incidence is smaller than in networks with no 
correlations. In large social networks this may lead to a 
difference of 15% to 20% of infected people for moderate 
values of A. On the other hand, we have found that the 
diseases are longest-lived in assortative networks. Addi- 
tionally, we stress that the method employed here can be 
used to solve other epidemic-like models in networks with 
any correlations such as the SIS and rumor spreading 
|23| models, provided that they can be described through 
mean- field rate equations. Finally, we point out that re- 
sults for disassortative networks will appear elsewhere. 

Y. M. thanks A. Vazquez for useful discussions. Y. M. 
acknowledges financial support from the Secretaria de 
Estado de Education y Universidades (Spain, SB2000- 
0357). This work has been partially supported by the 
Spanish DGICYT project BFM2002-01798. 



[1] 
[2] 

[3: 

K 

[5 
[6 

[7; 

[§ 
[9 
[10 

[11 

[12 
[13 

[14; 

[15 

[16 
[17 

[18] 



D. J. Watts and H. S. Strogatz, Nature 393, 440 (1998). 
S. H. Strogatz, Nature (London) 410, 268 (2001). 
S. N. Dorogovtsev and J. F. F. Mendes, Adv. Phys. 51, 
1079 (2002). 

A.-L. Barabasi, and R. Albert, Science 286, 509 (1999). 
R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 

(2002) . 

M. Faloutsos, P. Faloutsos, and C. Faloutsos, Comp. 
Com. Rev. 29, 251 (2000). 

R. Pastor-Satorras, A. Vazquez, and A. Vespignani, 
Phys. Rev. Lett. 87, 258701 (2001). 

R. Cohen, K. Erez, D. ben-Avraham, and S. Havlin, 
Phys. Rev. Lett. 85, 4626 (2000); 

D. S. Callaway, M. E. J. Newman, S. H. Strogatz, and 
D. J. Watts, Phys. Rev. Lett. 85, 5468 (2000). 
R. Albert, H. Jeong, and A.-L. Barabasi, Nature 406, 
378 (2000). 

R. Pastor-Satorras, and A. Vespignani, Phys. Rev. Lett. 
86, 3200 (2001). 

M. E. J. Newman, Phys. Rev. Lett. 89, 208701 (2002). 
P. L. Krapivsky and S. Redner, Phys. Rev. E 63, 066123 
(2001). 

A. Vazquez, and Y. Moreno, Phys. Rev. E 67, 015101(R) 

(2003) . 

M. Boguna, R. Pastor-Satorras and A. Vespignani, Phys. 
Rev. Lett. 90, 028701 (2003). 

M. E. J. Newman, Phys. Rev. E 67, 026126 (2003). 

Y. Moreno, R. Pastor-Satorras, and A. Vespignani, Eur. 

Phys. J. B 26, 521 (2002). 

M. E. J. Newman, Phys. Rev. E 66, 016128 (2002). 



[19] We assume that the conditions for the set of equations 

to be meaningful are satisfied |Tb1 |. 
[20] J. B. Gomez, Y. Moreno, and A. F. Pacheco, Phys. Rev. 

E 58, 1528 (1998). 
[21] Y. Moreno, A. M. Correig, J. B. Gomez, and A. F. 



Pacheco, J. Geophys. Res. B 106, 6609 (2001). 
[22] A. Vazquez, and M. Weigt, Phys. Rev. E 67, 027101 
(2003). 

[23] Y. Moreno, M. Nekovee, and A. Vespignani, in prepara- 
tion. 



